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Abstract. We present an approach to the dynamics of interacting particle systems, 
which allows to derive path integral formulas from purely stochastic considerations. 
We show that the resulting field theory is a dual version of the standard theory of 
Doi and Peliti. This clarify both the origin of the Cole-Hopf map between the two 
approaches and the occurence of imaginary noises in effective Langevin equations for 
reaction-diffusion systems. The advantage of our approach is that it focuses directly on 
the density field. We show some applications, in particular on the Zero Range Process, 
hydrodynamic limits and large deviation functional. 
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1. Introduction 

1.1. Generalities 

Many problems of current interest in statistical physics of out of equilibrium systems 
involve strongly interacting particles exhibiting non trivial collective phenomena. One 
example is that of supercooled liquids, where the dynamics slows down dramatically 
as the glass transition is approached, due to the increasingly collective nature of the 
dynamics, see e.g. pp. Another example is given by systems of diffusing particles that 
branch and/or annihilate; depending on the relative strength of these effects a variety 
of non-equilibrium transitions and anomalous scaling behaviour appear [2]. A third 
example is provided by systems driven out of equilibrium by external sources. A 
celebrated example is the one dimensional asymmetric exclusion process for which one 
finds phase transitions between different non equilibrium steady states. 

Developing theoretical techniques for such difficult problems is of great importance, 
in view of the diversity of situations in which they appear. 

A natural framework to study these collective phenomena is field theory which 
has been crucial in understanding equilibrium phase transitions. In the context of 
non-equilibrium systems, it has been already applied successfully to reaction-diffusion 
systems. It have been also crucial to get an handle of strong or intermediate coupling 
problems where no perturbative technique is at disposal. Two examples are the 
application of the Exact Renormalisation Group to the pair contact process [3J Hj and 
the Mode Coupling Theory of the glass transition [5]. 

A field theoretical formulation of interacting particle systems which has become 
standard is based on the Doi-Peliti formalism (DP)) [6j [7j. Starting from a second 
quantization representation of the Master equation, one obtains, after a rather elaborate 
coherent state representation, a field theory representation in terms of two fields </> and 
4> (see below). This has furnished the starting point of a very large number of studies 
[3], including Exact Renormalisation Group calculations. However, besides its intrinsic 
difficulty, the formalism is not transparently related to stochastic equations for the 
particle evolutions. Actually, the action of the field theory corresponds to a reasonable 
looking Langevin equation for the density of particles, except that the noise is often 
complex or even pure imaginary! This suggests that the field 0, despite its superficial 
resemblance with the density, in fact lacks a direct physical interpretation |8J. Other 
difficulties arise when one wants to treat systems of hard core particles or of particles 
with non-trivial diffusion constants. The relationship with stochastic equations on the 
particle trajectories is particularly important to study the hydrodynamic limit and make 
a connection with the large deviation functional techniques developed in [9j [101 (TTJ. 
The aim of this article is to discuss in full detail how stochastic questions on particle 
trajectories are related to field theory. We will unveil a dual version of the DP field 
theory that is naturally related to stochastic equations. This will shed new light on the 
underlying structure of DP field theory and we allow us to re-obtain some recent results 
from a different perspective, e.g. the large deviation functional of Bertini et al [11] We 
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will also present some new applications, e.g. we will derive the stochastic equations 
characterizing the dynamics of the Zero Range Process. 

1.2. Issues and questions through a simple example 

Let us illustrate what are the main questions and issues we want to address focusing 
on a simple example: particles A diffusing on a lattice with diffusion constant 7 and 
coalescing (A + A —>■ A) when they meet on the same site with rate A per unit of 
time. Following the DP formalism, which will be detailed below, the average density of 
particles in the systems is given by the average of a field in a path integral calculation 
with action 

S = f 1 dt f d D x (-d t <f> + 7A0) + A0 2 (1 - (f) 2 )} + boundary terms, (1) 
hi J L J 

after the extra field cf> has been integrated out. The quadratic part comes from diffusion, 
while the rest comes from the pair coalescence. After a shift (j) — > <fi+ 1, the above action 
becomes identical to the one obtained through the Martin-Siggia-Rose-DeDominicis- 
Janssen technique from the Langevin equation: 

d t <f> = 7A0 - 2A0 2 + 77 (r](x,t)r}{x',t')) = -2\<p 2 (2) 

The first two terms of the RHS are exactly what one would expect naively: a diffusion 
term plus an annihilation term. The problem is that the white noise has a negative 
variance implying that rj is purely imaginary! This problem, first observed in [12] 
demonstrates that the field is not equal to the density field, as can be seen from a direct 
computation. Furthermore, it suggests that either the physical sources of fluctuations 
in the system are not Gaussian or <fi is not a direct probe of fluctuations, or both. Still, 
stochastic equations for the density field certainly exist. Actually, in many cases one 
starts from a phenomenological stochastic equations to get the field theory and not the 
other way round. But then, understanding the relationship between stochastic equations 
for the particle trajectories and the DP field theory is crucial. This is the main aim 
of our article which is based on our joint recent work [5] with A. Andreanov and J.-P. 
Bouchaud. 

The main questions we want to address are: 

• What are the stochastic equations governing the evolution of particles? 

• How are they related to field theory? In particular, starting from these stochastic 
equations and using the MSRDJ technique what type of field theory is obtained? 
How is this field theory related to the standard DP one? 

• What is this field theory useful for? 

In the following we will answer the first two questions in detail. As for the third, 
we will show some applications, discuss some other promising ones and hope that the 
readers will find new ones. 
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1.3. Examples of systems of interacting particles 

In the following we will introduce different classes of systems to which our theoretical 
approach applies and on which we shall focus on in the next sections. 

1.3.1. Diffusion on a lattice and the Zero Range Process The simplest example is that 
of non (or very weakly) interacting particles diffusing on the hypercubic lattice A D . 
At each infinitesimal time step (t, t + dt) a particle jumps to a nearest neighbor with 
probability jdt/z where z is the site connectivity and 7 the diffusion coefficient which 
may depend on the position. 

The situation gets considerably more complex when 7 is made explicitly dependent 
on the number of particles at each site or in the neighbourhood. Such is the case in 
the Zero Range Process (ZRP), where a particle at site i jump to any neighbour with 
probability / ~fu(n i ), where is occupation number of site i and u(n) a function which 
vanishes for n = 0. This simple dynamic rule, although leading to factorized steady 
states, leads to interesting phenomena such as Bose-Einstein condensation [13]. Many 
variants exist, with several species [II], or diffusivity depending only on the target 
site [15]. 

1.3.2. Interacting particles Another class of interesting models is that of point particles 
interacting via some potential, which for simplicity we shall consider to be pairwise. 
This includes lattice systems where sites cannot be occupied simultaneously by several 
particles, provided this constraint is respected by the initial condition. It also includes 
systems like gases, simple liquids or crystals, when molecules can be approximately 
considered as pointlike. The Hamiltonian of the system reads: 

H = ^2v(xi-Xj). (3) 

i<j 

There are several possible ways of modelling equilibrium dynamics for such systems, 
which are equivalent in the continuum limit. We will focus in particular on the following 
two: 

(i) Particles jumping on lattice from sites to sites, using the Metropolis rule for jump 
acceptance, 

(ii) Particles in free space obeying the Langevin equation: 

-jjfc = ~Y. v ( x i ~ x i) +r li, (4) 

where rji is a Gaussian white noise with variance {rjiityVjit')) = $i,jfi{t ~ 0- 

If one chooses v(x) = for x > 2a and v(x) = 00 for x < 2a and choose initial 
conditions such that all particles are at least at distances 2a, the system is equivalent 
to hard spheres of radius a. Furthermore, if the particles evolve on a hypercubic lattice 
with spacing 2a, this is equivalent to condition of single site occupation. In the so 
called Exclusion Process, Symmetric (SEP) or ASymmetric (ASEP), particles diffuse 
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freely on a lattice with the restriction that they cannot overlap and (in higher than one 
dimension) must avoid each other in order to cross. The system exchanges particles 
with some infinite reservoirs at its boundaries, with some rates, or may be closed. As 
shown for instance in [9] for the one dimensional case, the role of boundary conditions 
is crucial. 

In these models, as in those described in the previous paragraph, the number of 
particles is conserved. In a continuous description, this is expressed by a continuity 
equation d t p = — V • J, where J is the local particle current. 

1.3.3. Reaction- diffusion systems In many situations, where particles may appear, 
disappear, or be changed into something else, the number of particles is not conserved 
and no continuity equation can be written, making hydrodynamic descriptions in 
principle more complicated. Paradigmatic models for such phenomena are given by 
particle assemblies with reaction and diffusion. For instance molecules deposed on a 
substrate, in a gas, or in a porous material, may diffuse and react chemically. The 
problem is in general to study the concentration of molecules in the non equilibrium 
steady state. Similar systems are those involving natural species, with birth and death 
processes, as well as predation, or epidemic spreading, where contamination occurs at 
contact. In general, such systems can be easily mapped onto systems with reaction- 
diffusion. If we denote by A and B two different species, several reaction processes are 
possible: 

(i) birth: -> A, 

(ii) death: A -> 0, 

(iii) coalescence: 2A —>■ A, or more generally mA —>■ nA, with m > n, 

(iv) contamination: A + B — > 2A, or pA + qB — > p'A + q'B, with p + q = p' + q' and 
q' < q, 

(v) transmutation: A — > B 

(vi) death at contact: A + B — > A, 

the list being not exhaustive. Generically, when they contribute, these processes occur 
at some rate A per particle per unit of time. 

A common simplification is to allow multiple occupation of sites and to consider only 
local reactions occurring at the same site. One can in principle consider particles with 
hard core repulsion or interactions occurring when the particles occupy neighbouring 
sites or compact clusters. 

2. Stochastic process and field theory for interacting particle systems 

In the following we show how to derive the field theory directly from the very definition 
of the stochastic process. This will bypass completely all the technical machinery based 
on coherent states used for DP field theory and it will make clear the relation between 
the stochastic process and field theory. 
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2.1. Basic processes 

2.1.1. Particle disintegration Let us start with the simple situation of a single site, 
occupied at initial time to = by iVo particles. The time is cut in intervals [t^ti+i] 
(z = 0, • • • , N— 1) of length dt, during which each particle disintegrates with probability 
A dt. Our goal is to characterize the dynamic distribution of the number of particles 
n{t) during the whole time interval [0, T]. To this aim, we introduce an auxiliary "jump 
process" J{U) = n(t i+ i) — n(ti) and a conjugated field n(ti), which will probe the 
generating function Z[{n(ti)}] = ( e dt ^ h{u)J{u) ). For a reason which will be clear below, 
we temporarily forget that J is related to the variations of the number of particles. Its 
statistics is very simple: J(U) = —1 with probability n(ti)Xdt and J{tj) = with 
probability 1 — n(tj)\dt. This gives 

Z[{n(U)}] = I] (l - n(U)X dt + n(U)\dt e ^ u) ) (5) 



when N — > oo and T = N dt is kept fixed. Following the standard approach developed 
by Martin, Siggia, Rose, Janssen and de Dominicis (MSRJD) [IB], the average of any 
observable 0[{n(tj)}], with fixed initial conditions is given by the following formula: 

(O) = ^(J {dn(U)} 0[{n(U)}] n^K^+i) -n{U) ~ J{U)))j (6) 

i 

where Z is such that (1) = 1. Here, the number of particles is ensured to be an integer, 
as well as its variations, thanks to the delta functions and to the initial conditions. 
Remark that no infinitesimals like dt appear explicitly yet, as they are hidden in the 
distribution of the J's. 

Using imaginary Fourier representations of the Dirac deltas, the average in ([6]) 
becomes 

(O) = ^(J {dn(ti)dn(ti)} 0[{n(ti)}} e^AWWM)^)-^))) (7) 

- / {dnitjdhiti)} 0[{n(ti)}} (eEi^H^O-n^J+^^Ei^)^** 5 - 1 ])). 
Z J 

Taking the continuous time limit (i.e. dt — > 0), one gets that the probability of observing 
a given path {n(t')}t'e[o,T\ is 

±J{dn( tl )}e~ s ^\ (8) 

with 

£[{71, n}] = - £ dt (-h(t)d t n(t) + Xn(t) [e _A W _ 1] ) . (9) 

The action ([9]) is made of two parts. The first one, with time derivative, is a "kinetic" 
term, which role is to fix the now infinitesimal increments of the field n. The second one 
comes from the jump process J and encodes the dynamic fluctuations. This structure 
in two parts is generic, and thus the steps made above, which consist in obtaining the 
probability distribution of the paths from the generating function Z[{n(tj)}] will be 
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skipped in the other examples. Remark here that the introduction of the field n allows 
to use real values for n instead of integers only. 

The generating function generates cumulants of the distributions of the variation 
of particle numbers. In particular, the formula 

5Z[{n}} 



(d t n(t)} = 



(10) 

n(t)=0 



8h(t) 

will be useful for the derivation of hydrodynamic equations. 

2.1.2. Particle creation Another simple process has been described briefly in the 
introduction. The calculation follows the lines of the previous paragraph, and gives 

S[{n,n}} = dt (n(t)d t n(t) + A [e" w - l]) . (11) 

2.1.3. Diffusion The diffusion process can be either seen as the motion of a particle 
performing a sequence of jumps on connected sites, or as an exchange of particles 
between nearest neighbours. As we are interested in the density fluctuations, we adopt 
the second point of view, which amounts at considering only two neighbouring sites 
1 and 2, with initial occupations ni(0) and n 2 (0). Particles may hop back and forth 
with rate '-/dt. The (integer) numbers of particles on the two sites at time t are n\(t) 
and n2(t). The variation of n& between t and t + dt will again be noted dJk{t), while 
particles hop from 1 to 2 with rate Wyi and from 2 to 1 with rate W^i- Of course, 
dJ\(t) and dJ2(t) are strongly correlated since a particle leaving site 1 lands on site 2 
and vice- versa. More precisely, dJi(t) = —dJ 2 (t) = +1 with probability n2(t)W2\dt, 
dJ\(i) = — dJ 2 (t) = — 1 with probability ni(t)Wi 2 dt, and dJi(t) = c?J 2 (t) = otherwise. 
As before, the dynamical action is obtained through use of the generating function, 
leading to 

S[{n,h}\ = - j dt{-h 1 d t n l -h 2 d t n2 + n l W l 2{e h ^ - 1) + n 2 W 2 i(e fll - fl2 - 1)}. (12) 
On a lattice, the total MSRJD action reads: 

S[{n,n}] = - f dt l-^nAni + ^n^e"^ - 1) i , (13) 

where the brackets restrict the summation to pairs of nearest neighbours. Later, we 
shall use this action to get continuum limits. 

2.2. Interacting particles and Dean equations 

Let us consider particles diffusion on a lattice with fixed hopping rate, but feeling 
a potential Ui(t) at site i. This potential can be external and time dependent or 
can also be due to the interaction of particles at site i with all others, in particular 
Ui(t) = a D v(a\i — j\)ni(t)rij(t). A natural stochastic dynamics is to accept hops 
to neighbouring sites with the heat bath rule, that is by choosing 
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where (3 = I/7 is the inverse temperature. In the continuum limit Wij can be expanded 
in powers of a(i — j) and Langevin dynamics is recovered. 

Before doing so, it is important to define the continuum limit properly. We start 
with off-lattice particles in a box of very large size L with given boundary conditions. 
We then devide this box into many tiny boxes of size a at position x = ai. We use the 
number ni of particles in this box to define the density p x = riijaP and its conjugate 
p x = hi. We choose a very small in order to expand in powers of a. At the lowest useful 
orders: 

gftf-fi, _ 1 = gOe^Vfii+^eyV)** _ x + ^2) ^ 



aeij ■ Vhi + c) 



a 2 

(ey ■ + (ey • Vfiif + o(a 2 ), (16) 



where is the unit vector pointing from i to j. One can also expand the expression in 
and get 

W^ = 7--aey-VC/i + o(o). (17) 

In addition, using x = ia, one can replace aP J2i by / d D x. Terms of order a vanish 
by symmetry, and thus keeping terms of order a 2 and rescaling t by a 2 , one gets the 
following action: 

S[{n, h}] = - J dtj d D x {p x {-d tPx + ^Ap x + V ■ [ Px VU x ]) + 1Px {Vp x f} (18) 

This action can also be obtained from a MSRJD treatment of the following Langevin 
equation: 

d t p x = -V • J x , (19) 
with a fluctuating current 

J x = -7V p x - p x VU x + v/p^, (20) 
where £ x (t) is a Gaussian field with variance 

(Ut%>(t')) = 2 1 6(x-x>)6(t-t'). 

The deterministic part of the current can be expressed in the form — AzV|^-, with a 
free-energy functional 

F[{P}\ = J d D x (Tp x ln(a D p x ) + p x VU x ) . (21) 

This stochastic differential equation (SDE) for the density field was first derived by 
Dean [17] directly from the Langevin equation for interacting particles by using Ito 
formula. The free energy T is the same as that from mean-field theory. Our derivation 
gives further evidence that no coarse-graining is needed in order to get such free energy. 
In fact the interaction U used here is often obtained as an effective interaction (eg. 
for colloidal systems) and thus the free energy functional is actually different from the 
mean-field one, which contains the microscopic potential instead. 
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2.3. Newtonian dynamics 

Now let us consider systems of particles with Newtonian dynamics and interacting via 
a pairwise potential v. The equations of motion read 

dtXi = V - (22) 
m 

dtPi = -J2^ v ( x 'i~ x j) ( 23 ) 

We introduce the density field in phase space: 

p(x,p,t) =^2S(xi{t)-x)S(pi(t)-p). (24) 

i 

It is easy to show that it obeys the deterministic equation 

d t p(x,p,t) = ——d x p(x,p,t) + / d D x'd D p' p(x',p', t)W(x — x') ■ d p p(x,p, t), (25) 
m J 

all randomness being in the initial conditions. One can write path integral formulas as 
before, introducing p(x,p, t) fields. If one applies the reverse map of fl60l to these fields, 
one gets the formulas originally obtained by Doi [6]. 



2-4- Exclusion processes 

Dealing with hard cores is in general difficult. A first attempt was done in [18] Particles 
with hard core repulsion may be modelled in two ways. The first one is to add a pairwise 
interaction which is infinite if the distance between the particles is less than their radius 
and vanishes else. The second one is to include a constraint in the hopping rate of the 
particles on the lattice. From the stochastic approach which has been developed here, 
the later is easier to deal with. We thus start with particles on a lattice with initial 
occupation number of site i rij(O) G {0; 1}. The exclusion constraint is propagated by 
the dynamics, if the hopping rate from i to j is of the form 

UiWij = AijTii(l - rij), (26) 

which vanishes if site % is empty or site j is occupied. For sake of clarity, we choose 
Aij = 7, which corresponds to symmetric exclusion processes. Inserting (1261) into (fT3]) , 
expanding like in 015|) and using nj = + ae^ ■ Vrij + o(a), the action describing the 
dynamics becomes in the limit a — ■> 0: 

S[{n, h}) = -gT d J dt J d D x {-n x d t n x + jn x An x + jn x (l - n x ) (V^) 2 } , (27) 

after rescaling of the microscopic time t — ► t/a 2 . This looks like the hydrodynamic 
equation obtained rigourously by Bertini et al. [11] for symmetric exclusion processes 
- althought the time rescaling is different. However, our approach suggests that this 
equation is valid beyond the hydrodynamic limit. In addition, although we started from 
n x G {0; 1}, the measure gives also non-zero weights to paths such that n x (l — n x ) ^ 0. 
As we shall see, this is related to the fact that the hydrodynamic limit corresponds to a 
saddle point of the dynamic action which is located on continuous paths, and that the 
dynamics on slightly smaller scales is given by fluctuations around this saddle point. 
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2.5. Initial conditions 

Up to now, we have not specified the initial conditions and considered only the dynamical 
changes of the density. In many cases, initial conditions are forgotten after some 
transient time and a stationary state is obtained. However, in the case of exclusion 
processes, our stochastic approach to the dynamics makes it mandatory - at least 
for consistence - to specify initial conditions with the good occupation numbers, in 
order to insure that the paths supporting the measure are consistent with the dynamic 
rules. We assume a factorized initial state with particles on site i, drawn with a 
distribution Pj(<&). It is useful to define the generating function gi(z) = In (j2 q P i (q)z q ^j. 
In order to specify the initial conditions, we add a term (n* 5 (qi — rij(0))). Using integral 
representations of the <5's, this gives an extra term to the dynamic action, after averaging: 

5i[{n(0),n(0)}] = £n,(0)n,(0) - J> (e*<°>) . (28) 

i i 

A quite generic case is that of Poissonian initial state, with Pi(n) = -^fe _Po ' i , which 
leads to 

5![{n(0),n(0)}] = £ n,(0)n,(0) (e*<°> - l) . (29) 

i i 

If hard cores prevent several particles to sit on the same site, one may choose independent 
Bernoulli distributions with average po t i, leading to 

Si[{n(0), n(0)}] = ]T ni (0)n t (0) - ]T In (l - p 0ji + p ,,e^ (0) ) . (30) 

i i 

2.6. Boundary conditions 

It may happen that changing boundary conditions changes drastically the dynamical 
behaviour. One dimensional exclusion processes are popular examples of such 
phenomenon. Periodic boundary conditions are taken into account by imposing space 
periodicity on the fields. Taking into account open boundary conditions is more model 
specific, and we explain how to do it on the example of an exclusion process in one 
dimension. We assume that the system is coupled to two reservoirs at its ends. Particles 
are added at boundary points and N with rates a and 7, and removed with rates (3 
and 5. As the number of particles is preserved inside the system, its fluctuations are 
governed by exchanges of particles with the reservoirs. This gives an extra term to the 
dynamical action: 

S B [{n ,n N ,h ,n N }] = - J ' dt {a(l-n (t)) (e"°^ - l) (31) 
+ 7 (1 - n N (t)) _ 1) + ^(t) ( e -M*) _ 1) + S n N (t) (e""^ - l)} . 

3. Mapping between Doi-Peliti and stochastic approaches 

This section is devoted to show mapping between the above approach and the more 
standard one based on the DP method. For self-consistency, we start by recall briefly 
the derivation of path integral formula using the latter method. 
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3.1. Doi-Peliti path integral formulation 

Let us consider particles with annihilation A + A — > on a single site at rate 2A, which 
will provide a benchmark to show our method. The initial number of particles is hq. 
The probability of having n particles at time t evolves according to the master equation 

d t P(n, t) = X (n + 2)(n + l)P{n + 2, t) - X n{n - l)P(n, t). (32) 

We map this onto a imaginary time Schrodinger equation by introducing a Fock space 
generated by n-particle states \n) and introducing the state ket of the system at time t: 

m))=Y, P (n,t)\n). (33) 

n 

We also use the lowering/ raising operators a and a + according to the standard definition 

a|0) = (34) 
a + \n) = \n+l) (35) 
a\n) = n\n — 1), (36) 

which verify the commutation relation [a, a + ] = 1. The master equation ( |32l) now 
becomes 

d t \m) = -H\m), (37) 

with 

H = X (a +2 - l) a 2 . (38) 
We obtain the state ket at any time: 

Mt)>=e-**M0)>, (39) 

and 

P(n,t)= K (40) 
n! 

The ra-particles states verify (n\m) = 8 n , m n\. As a consequence one can define coherent 
states 

' V^ a+ |0) (41) 



/7T 

and 



1 e-^(0\e^ a . (42) 

7T 



The coherent states thus define the (over) complete relation: 

1 = J d(f)d(f)*\(j)}((j)\. (43) 

The average of an observable A(n) is given by 

(A) t =J2 A (n)P(n,t) (44) 

n 

= (P\Mm) (45) 
= (P\Ae-™mO)), (46) 
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where 

i = EW#H (47) 

n n - 

and the "projection state" (P\ = (0\e a is a left eigenvector of a + with eigenvalue 1. 
Using the Trotter formula 

M 



e~ Ht 



(l - Hdt + o(dt)j , (4? 



with M = t/dt, and inserting ( 1431) in between all factors, one obtains a product of terms 
of the form 



It is easy to show that 

((pk+i\H\4>k) = (<f>k\<f>k)H(<t>i + i,<f>k), (50) 

where if (</>*, </>) is obtained by replacing respectively a + and a by 0* and in the 
expression of : H : in terms of a + and a. The expression : H : is obtained from H 
by normal ordering, i.e. commuting all operators until all a's are on the right. In this 
example (|38|) has already been normal ordered, and thus 

#(<^,0) = A(0* 2 -l)0 2 . (51) 

We thus get: 

. M 

(A)t = J II d ^d<Pl (P|i|0A/)(0o|^(O)) (52) 



H\<f> k )' 



x 



k=0 

' M \ 

n<0fci0fc-i)(i-d^(^>fc-i)) ■ 



The operator A can always be written in the normal ordered form A(a + , a) =: A(a + , a) :. 
In addition, using the identity [e a , f(a + )} = f(a + + l)e a , we get 

(P\A\<l> M ) = A(l,<f> M )e+". (53) 

In addition: 

<0oh/KO)>= T,P(q)(M(a + ) q \0) (54) 



= EHlW (55) 

= e 9( ^ } . (56) 
In the limit dt — > 0, we introduce the continuous time s = kdt and thus 

(A) t = J {d<P*}{d<P}A(l,<p t )e Ht)+9ir ° )+ fo dsH{r{s) ^ (57) 
Diffusion processes can be treated in the same way. The resulting hamiltonian is 

#dif [{</>*, 0}] = 7 E 0* " - h) (58) 
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Correlation functions at different times can be computed in the same way. For instance: 
(7i(ti)---n(tp)> = (P\a + af[{e- {t ^ tk - l)H(a ' a+) a + a}^2P(q)\q)- (59) 

k=2 q 

The average of the 0* field is 1, and in general the shift = 0* + 1 is carried out, in 
order to deal with a field with vanishing average. This also cancels the isolated term 
0(t) in dSTJ). 

3.2. The Cole-Hopf transformation 

We now show how the action obtained from purely stochastic considerations can be 
exactly mapped onto the one derived in the previous paragraph. The particle number 
operator is n = a + a. Thus, it is tempting to express the density field as p = 4>*(f). 
However, in order to go from operators a,a + to fields 0,0*, one must normal order. 
For instance : h 2 := (a + ) 2 a 2 + a + a, and thus in the field theoretic formalism, one 
should rather choose p 2 = (0*0) 2 + 0*0 instead of p 2 = (0*0) 2 . Before dealing with 
this ambiguity, we will proceed to a naive mapping between fields in the different 
formulations, and apply it to situations where it is correct. Let us define 

p i {s)=<j>*<j) i {s), (60) 
p t {s)= ln0*( S ). (61) 

The mapping ( 1601) has unit determinant and thus does not change the measure. In 
addition, it is valid both on lattice and in the continuum. Plugging (|60|) into (!58|) . one 

gets 

7E fdspi(s)(e^-^-i). (62) 

This is precisely the terms resulting from the calculation of the generating function in 
the stochastic approach. The extra kinetic term in ( 1571) becomes under pOj) : 

-V/ ds pi(s)d s pi(s). (63) 
: Jo 

This shows that the field defined as p = 0*0 is a "good" density field, as it formally 
corresponds to the one of the stochastic method. The nature of density fluctuations is 
thus hidden in the DP formalism. Thanks to the Cole-Hops mapping detailed above 
one obtains a dual version of the field theory which turns out to be the DP one as the 
reader can easily check for all processes considered in this paper. There is only one 
subtle point needed for this derivation that we shall discuss below. For concreteness let 
us come back to the reaction A + A — > 0. The non-kinetic bulk part of the dynamic 
action reads in the DP form: 

S a+ a,dp = - J*ds\<f> 2 (s) (l - 0*(s) 2 ) . (64) 

The equivalent term in the stochastic approach is: 

Sa+a,msrjd = ~ f ds Xp(s) (p(s) - 1) (e- 2 ^ - l) . (65) 
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However, if we use the mapping (1501) into (IM1) . one gets 

ds \p{sf (e- 2f,{s) - l) , (66) 

which differs from ( 1651) . In order to understand the origin of this difference, we recall 
that the fields and 0* are quantum fields, as a consequence the products of fields at 
the same time is a quantity that depends crucially on the underlying discretization. For 
example (p*(t + )(p{t~) ^ (p{t + )(f)*{t~) as it can be readily checked. In order to understand 
which discretization allows one to map the field theory derived in the previous section 
to the DP one we carry a mapping equivalent to (160)) at the operator level: 

a = e~ p+ p (67) 
a + = e p+ . (68) 

Operators p and p + have canonical commutation rules: 

[p, p + ] = a + [a, lna + ] 

= a + d In a + /da + (69) 
= 1. 

It is thus natural to use them as lowering/raising operators. We define a new set of 
Fock space generating vectors as follows: 

|0> = |0>, (70) 
|n) = (p+r|0>, (71) 
and corresponding left vectors: 

(0| =(P\ (72) 
(m\ = (0\p m . (73) 

We point out that p and p + are not hermitian conjugates, and thus neither are (h\ 
and We can redo the DP derivation of the path integral formulas, but using these 
operators instead of operators a and a + , writing 

H(a + ,a) = H stoc (p + ,p). (74) 

Normal ordering with respect to p + and p and replacing respectively p + and p with p 
and p, one gets the equivalent of (|57|) . where the term at time s is 

- p(s)d s p(s) + H stoc (p(s),p(s)). (75) 

Initial conditions are also taken into account: 

£P(g)</«= £P(g)e«o* 

= e a ^\ (76) 

Gathering terms from fl75]) and (IT5"|) . we get exactly the dynamic action derived from the 
purely stochastic approach. From the point of view of the time discretization the above 
derivation clarify what was the origin of the discrepancy between (I65p and (IB"6"|) . The field 
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theory derived directly from the stochastic equation gives, after fields transformation, 
p{t)p{t) = lim £ _ tO 0*(£ + 4e)0(t + 3e)0*(t + 2e)0(t + e). Instead within the DP field theory 
the corresponding term is 0*(t)0*(t)0(t)0(t) = lim e ^ 0*(t+4e)0*(t + 3e)0(t+2e)0(t + e). 
As a consequence, in order to map one onto the other one has to use also the same time 
discretization. 

We remark that we have not tried to justify the existence of the operators p and 
p + . The change fl67|) may not be well defined at the operator level. However, it is a 
posteriori justified by the fact that it gives exactly the action obtained from another, 
stochastic, approach. Furthermore, the same result can be obtained alternatively just 
by discretizing the field theory. 

As a consequence, although path integral formulas can be obtained both using the 
simple MSRJD stochastic method or the more abstract DP one, the previous discussion 
highlights that care must be taken in the interpretation of the field theory. Indeed, 
depending on the physical situation, the sets of fields 0, 0* and p, p are more adapted 
than the other. However, different sets are adapted to different approaches and going 
from one to the other using fl60|) may be dangerous. In fact this mapping has been often 
referred to [3[ [121 [13 EE ED], in a rather loose way, in particular without attention to 
the fact that the commutation rules of the fields are important. 

3.3. Physical origin of fluctuations 

Let us discuss in more detail the origin of the difference between the two field theories 
obtained in the previous sections. The noise generated by the stochastic dynamics is 
intrinsically Poissonian at the microscopic level. For example, consider particles created 
on a site (empty at t — 0) with rate A per unit of time. The generating function of N(t), 
the number of particle at time t is f(s) = (e sN ^'). It is a straightforward exercice to get 



in the limit dt — > keeping t fixed. Similar expressions, are obtained for diffusions 
and more complicated processes. The intrinsic Poissonian nature of the noise, i.e. 
the fact that with very small probability variables change of a finite amount, makes 
the logarithm of the generating function a complicated function of the argument. In 
particular it leads generically to exponential terms. This is the main reason why the 
action of the field theory in the p, p variable is complicated even for simple processes. 
On the other hand, the field theory in the 0, 0* variables is simple for simple processes as 
creation, annihilation or diffusion. In particular for diffusion the action is quadratic. The 
drawback is that is not the density field. This is the reason why stochastic equations on 
doe not make sense physically and lead to inconsistencies, as the presence of imaginary 
noises. 

The usefulness of each of the two field theories depends on the particular physical 
problem one is interested in. One focuses directly on the physical variables but it 
is already complicated (non-Gaussian) for simple processes. The DP one is instead 
simple in these cases; however computing physical observables may be quite complicated 




(77) 
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because they correspond to high order multipoint correlation functions. Thus, for 
renormalization group treatment of simple processes the DP one seems and has been 
proven to be very useful. On the other hand, in cases where one is interested in 
two point density correlation functions the other field theory seems to be more useful 
because it allows to develop approximations directly for these functions avoiding to deal 
with high order multipoint (more than two points) correlation functions which is in 
general a quite cumbersome task. Finally, we remark that on sufficiently large length 
scales the noise becomes Gaussian. In this case the field theory on p, p leads directly 
to hydrodynamic descriptions away from critical points. Hydrodynamic descriptions 
correspond to coarse-graining the dynamics on mesoscopic scales much bigger than the 
correlation length, with corresponding rescaling of time. From coarse-graining results a 
law of large numbers under the form of a large deviation theory, providing the probability 
of occurence of coarse-grained trajectories, or of density profiles in the steady state. 

4. Applications 

We now give some applications of the formalism which we have explained in details. We 
start with the zero range process for which we derive continuum stochastic equations. 

4-1. Zero Range Process 

Zero Range Processes (ZRP) provide a simple model for the study of particles undergoing 
an out of equilibrium process. Indeed, the dynamic rules, although quite simple may 
lead to a rather reach phenomenology. In addition, in the simplest cases, the steady 
state measure can be computed and is factorized. For a detailed review, see [T3] . 

4-1.1. Single specie ZRP The model is defined as follows. Particles hop on a lattice 
with a rate u(rii(t)) depending only on the number of particles at the considered site i 
at the moment of the jump ; u is a function which vanishes at (no motion from site 
i if no particle). Here the continuum limit is obtained in a way slightly different from 
before. The lattice spacing is a, and is sent to zero, while the density is defined as before 
from the particle number. It is now a straightforward exercice to generalize what has 
been done previously to get the following action: 



with u a (p) = a D u(a D p), after rescaling of the microscopic time t — > t/a 2 . Integrating 
out the p field like previously, one gets that the dynamics of the density field follows the 



S[{p,p}] = - JdtJ 





SDE 



d t p. 

J x 



V ■ J, 



X 



(79) 
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where £ has the same statistics as in (|2"0|) . The steady state free energy can be 
obtained from the Fokker-Planck equation. We introduce the instantaneous probability 
distribution of paths: 

n{p},t] = (uHp x -p x (t)yj, (si) 

which obeys the following Fokker-Planck equation: 

dtV[{p},t] = - (J^,Au a (p x )V[{p},t] + V ■ (u a (p x )VV[{p},t})^ , (82) 
Where we used the notation (f,g) = J d D x f x g x . Writing formally 

Au a (p x ) = V • ^(fe)vSJ , (83) 
Eq. ([S2D reads 

^ = - ( J-> v • ( M ^) v ^^) v + v • K(p-) v ^)) • 

It is straightforward to check that the steady state distribution is obtained from the 
solution T of (1531 . which is 

T\{~p\\ = I d D x [ PX dr\nu a (r), (85) 

J Jl/a D 

via 

P[{p}]ocexp(-^[{p}]) (86) 
This is the continuous analog of 



J4) 



i 



where u\ n = Ilfc=i Uk- The possibility to solve (!83|) by direct integration is related to 
the fact that the steady state measure is factorized. In general, this measure is not 
factorized, and it is more difficult to read T in the Fokker-Planck equation. We will 
illustrate such difficulty in the context of the two-species ZRP. 

We want to end this section with a remark on the continuum stochastic equations 
derived for the ZRP. They are valid on scales / such that a « I. On the other hand there 
is no assumption on the correlation length £. Coarse-graining the previous stochastic 
equations on scales much larger than £ would lead to hydrodynamic equations (see the 
next sections). The main interest of eqs. (1791) is that they provide at the same time 
a continuum description and can describe critical and out of equilibrium properties of 
ZRP, ie. they have the same status as the model A stochastic equations for the (non 
conserved) dynamics of the Ising model. 

Remark that we have taken spatial derivatives of the fields without justification of 
their smoothness. In order to make this more rigorous, one would have to consider the 
density fields as distributions acting on the space of smooth functions with compact 
support on the discrete lattice. This is in fact a rigorous but much less clear way of 
doing the coarse-graining which we refer to. 
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4-1-2. Two species ZRP This version of the ZRP involves two species A and B, with 
occupation numbers nf(t) and n B (t). The difference with the single-species case, is 
that the hopping rate from site i of particles of a given species is a function of the local 
number of particles of the other species, which we denote by u A {nf) and u B (n B ). In 
the continuum limit, the density fields verify the SDE 



dtpt = Artf (p£ , p B ) + V • \Jui(p*, P B X W X 

d tP B x = Au B ( P t p B ) + V • (vV(p^P B )£f) ■ (89) 
If the steady state is factorized, the following conditions must be verified: 

hxui(pt Px) = J~X + constant, (90) 
8T 



In u B {pt Px) = + constant. (91) 
This gives a necessary condition for factorization of the steady state: 

u B Mrp B J uAa{ ^ p f x) = uXptP B J UBa f p f x) - 02) 

This is the continuous analog of the condition given by Evans and Hanney for the model 
on-lattice 



4-2. Hydrodynamic limits 

We now show how hydrodynamic limits can be obtained within our field theory 
approach. Taking hydrodynamic limits consists in splitting the volume of the system 
into boxes equal sizes, with volumes a fraction of the total volume of the system. 
This amounts to using a coarse-graining length I much smaller than the system size 
L, but much larger than the lattice spacing a. In one dimension, boxes are thus located 
between xL and xL + I, with x G {0, l/L, 21 /L, ■ ■ ■ , (L — l)/L}. In the hydrodynamic 
limit, where both / and L go to infinity, x becomes a real coordinate in [0, 1] (this 
corresponds to dx = l/L). In order to avoid artificial showness, the time is rescaled: 
t — > t/L 2 . The hydrodynamic limit is in principle different from the much less controlled 
continuum limit taken before. There are in fact two possibilities. The first one is that 
the hydrodynamic limit exists, in which case the limits are connected ; this happens for 
instance away from criticality, where correlation lengths are finite. In this case the large 
time and length scale behavior of the continuum field theory leads to the hydrodynamic 
limit. The other possibility is that the hydrodynamic limit does not exist, while the 
continuum limit of the field theory still exists, although great care must be taken; 
this typically happens at criticality where the field theory in the continuum allows for 
renormalization group calculations. 

There are now standard routes for rigorously deriving hydrodynamic limits for 
particle systems on lattice, such as ZRP or exclusion processes. We will not repeat 
them, as they can be found in textbooks and recent papers. Instead, we will show 
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how our field theoretic approach naturally leads to hydrodynamic equations and large 
deviation functionals. Although the following results are not new, we think that the 
derivation we present gives an interesting new perspective on these problems. 

The basic idea underlying hydrodynamic equations is the existence of a law of large 
numbers for density profiles. Let us imagine a system with a finite correlation length 
£, coarse-grained using a coarse-graining length I ^> £. Inside all boxes, the density 
field is the sum of a large number of independent identically distributed processes and 
is thus Gaussian. It can thus be characterized by first and second moments, and a large 
deviation functional is easily deriveed in terms of the macroscopic transport coefficients 
- see last section. 



4-2.1. Zero Range and Exclusion Processes Let's start by a simple example, single 
species ZRP in one dimension. In such models, the stationary measure is factorized, 
which makes it possible and simple to compute local averages from the marginals of 
the total joined distribution. Coarse-graining on the scale / is equivalent to averaging 
out all fields while fixing the value of the density inside each boxe at each time to its 
average value ~p(x,t). This amounts to using constant conjugated field n inside each 
box. However, instead of restricting the path integral to piecewise constant ft fields, it 
is simpler to restrict it to slow varying h. More precisely, we restrict the measure to 
fields which significant variations occur on scales larger than I. The dynamical action 
reads when the lattice spacing a — l/L goes to zero: 



J dt {-UidtUi + a 2 u(n,j) {Vhi) 2 + Arc*] } 



(93) 



Now comes the fundamental hypothesis of local equilibrium. The variation of the total 
number of particle inside large boxes is only due to particle flows at its boundaries. If 
the size is large enough, one may assume that the time scale for any significant variation 
is also large and the system is locally at equilibrium. The local average density evolves 
slowly following the hydrodynamic equation which we want to derive, and we make the 
following approximation: 



y £ u{m{t)) [(Vn,(t)) 2 + An l (t) 

ieB x 

where 



y £«(«)*(*) [(Vn l (t)) 2 + An l (t) 

ieB x 



«(t) = «t)) 



(94) 



(95) 



i€B x 



is the local mean value of u(rii(t)). The assumption of local equilibrium also gives this 
mean value: 



= ( u ( n ))p x {t), 

where (-)p stands for the average using the marginal with mean p. 
It can be obtained from the grand canonical partition function 



(96) 



(97) 
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with p n = — {— : 

Inverting f[9"gj) defines R(p) = z. Then we get: 



u(n)- = z = R(p). (99) 
In addition, n being a very slow varying field, one has in one dimension: 

£ [(V^) 2 + Ahi] w / [(V^) 2 + Ap x ] , (100) 

where hi ~ p x inside the box of size / around x and V x is now the gradient at the 
intermediate scale x, i.e. the variation between two points at distance I divided by 
/. Carrying the same analysis with the term —J2i n idtni, and rescaling the time by 
a 2 = 1/L 2 the action finally reads: 

Shydro = -L J dt £ dx {p x d t p x + p x AR(p x ) + R(p x ) (y Px ) 2 } . (101) 

The result (11 Oil) is a good illustration of the emergence of a law of large numbers at 
the hydrodynamic level, where the empirical current has mean V-R(p) and mean square 
R(p). This can be seen either from a direct interpretation in terms of an effective 
Langevin equation or by using ffTU]) and second derivatives of the generating function. 
In particular ({TO]) gives the hydrodynamic equation: 

d t p = AR(p). (102) 

The action (II 011) gives the large deviation functional for density fluctuations, by 
integration on the field p. If one formally defines a(p) = —V ■ (i?(p)V), then the 
weight of any - coarse-grained - trajectory is given by 

P[{p}} « J {dp x (t)}e~ Sb y^ l{ ^ }] , (103) 

which gives 

P[{p}} oc e -^° m] , (104) 

with 

S hydm [{p}] =Ljdt J\d t p x - AR(p x ), a(j> x )~ l (d t p x - AR(p x ))). (105) 

Due to the L factor in front of the integral, the path integral is dominated by the 
hydrodynamic trajectory, which minimizes the large deviation functional. 

The case of ZRP is relatively simple as the stationary distribution is factorized, and 
thus u(n) is easily computed. In general, the whole product measure must be used. 

Starting from (127|) . the fluctuating hydrodynamic equation is obtained for exclusion 
processes, where a drift velocity v can be added: 

d t p = -V ■ (vp(l - p)) + Ap + V ■ (v/2p(l - p)/L 77) , (106) 

where 77 is as usual a normal Gaussian noise field. This equations contains initial 
symmetries of the microscopic model, i.e. particle-hole symmetry. 
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4-3. Systems with reaction processes 

Reaction processes can provide an example of hydro dynamic limits in which the time 
has not to be rescaled because there are no conserved quantities. 

Again, we take for canonical example the pair annihilation A + A — > 0, but without 
diffusion. Using the notations used until here, one gets 

A£ J dsmWtmis) - 1) (e- 2 "^ - l) « LXj dt£d D x(n(n - l))^ (s) 

x (e- 2Ms) -l) (107) 

and thus 

Sdyn ^-Ljds jf* d D x {\(n(n - l))^ (s) (V 2 ^ - l) - p x (s)d s p x (s)} . (108) 

The difficulty is still to compute the static average. For the simple example given 
here, there is no steady state except the empty one, but occupation numbers remain 
Poissonian at all times, and thus 

(n(n - = p(s) 2 . (109) 

As L is large, the measure is dominated by saddle points of the action. Thus the 
hydrodynamic equation is 

d t p x (t) = -2\p x (t) 2 e- 2Mt \ (HO) 

where p x {t) is solution of 

d t p x {t) = -2Xp x (t) (e- 2 ^) - l) . (Ill) 

This equation describes the behavior of the system on large length scales but timescales 
of order of one. There is no need of rescaling the time because the density is not 
conserved so the coarse-grained density still evolves on timescales of order one. 



4-4- Large Deviation Functional 

Recently there have been a lot of interest and works on large deviation functional for 
out of equilibrium driven systems. In the following we would like to show how this 
large deviation functional appears naturally within our framework. Our approach is not 
rigorous compared to the previous ones [TH [9] . On the other hand it shows clearly in our 
opinion the key ingredients and it can be easily generalized to systems more complicated 
than the ones considered up to now. In the following, for the sake of simplicity, we will 
focus on one dimensional driven stochastic lattice gases. The driving can be due either 
to the boundary conditions or to a force not deriving from a potential. 

We shall focus on the probability that the system follows a given path {n(x,t)} 
in configuration space. This can be obtained formally by integrating out the auxiliary 
field n but this is in principle not feasible. However, if one is only interested in the 
probability of a given path after coarse-graining then the tasks simplifies a lot. Let us 
rescale time and length scales in the hydrodynamic way described before and focus on 
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the probability of smooth paths. This can be obtained with a construction very similar 
to the one used to obtain Feynman path integral in quantum mechanics: 

P[{p(x,t)}]~ n PA X [p;t^t + At] (112) 

Ax, At 

where P& x [p;t — > t + At] is the probability inside a small coarse-grained box of the 
evolution of the density on a small coarse-grained time (p is the notation for the coarse- 
grained density). If the hydrodynamic limit exists then each Pa x [p', t — > t + At] can be 
replaced by its hydrodynamic expression because although the global fluctuation of the 
density can be large, on each small (coarse-grained) box they are small. This expression 
equals the corresponding current distribution P& x [J;t t + At], up to a constant 
Jacobian. Repeating the procedure explained previously for the hydrodynamic limit 
one obtains a field theory such that the auxiliary field can be integrated out leading to: 

P Ax [J ; t^t + At]~ J Vpexp[-LAxAt(J - JM) 2 /x(p)] (113) 

where J av (p) = —D(p)^7p + x{p)E ls the empirical current; D,x are the transport 
coefficients characterizing the hydrodynamic limit, i.e. diffusion coefficient and mobility 
(L is the size of the system) and E is the external driving field. Putting together all the 
probability on small boxes one obtains the large deviation functional: 

P[{p(x,t)}] ~ J Vpexp [-L J dt(d tP + V • JMMp) ^tP + V • J av (/>))>] , (114) 

where n(p) is the inverse of the operator —V • (x(p)V). Note that this field theory 
corresponds to a stochastic equations for the density field: 

d t p(x, t) = -V • J av (p(x, t)) + V • [y/x(p(x,t))rj(x, t)] (115) 

where rj is a white noise in space and time, (r)(x, t)r)(x', t')) = 25 (t — t')5(x — x')/L, and 
the multiplicative noise must be interpreted in the Ito sense. 

This derivation although certainly not rigorous has the virtue of showing in a simple 
manner why only hydrodynamic transport coefficients matter for the derivation of the 
large deviation functional and why this is related to a stochastic equation which leads 
to a non-linear fluctuating hydrodynamics. It can be straightforwardly generalized to 
systems with a more complicated hydrodynamics, e.g. real fluids where density, energy 
and momentum are conserved. 



5. Conclusion 



The aim of this paper was to discuss and present a field theoretical approach to 
interacting particle systems. We wanted to show a dual version of the Doi-Peliti field 
theory that can be obtained from the stochastic process and that is directly related to 
stochastic equations. The advantage of this field theory is that it focuses on the physical 
density fields as we have shown in some applications, e.g. ZRP continuum stochastic 
equations, large deviation functional, hydrodynamic limits. We think that this approach 
will help to tackle difficult problems in interacting particle systems. 
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